Numerical Analysis of the Cahn-hilliard Equation and Approximation for the Hele-shaw Problem, Part I: Error Analysis under Minimum Regularities

In this first part of a series, we propose and analyze, under minimum regularity assumptions, a semi-discrete (in time) scheme and a fully discrete mixed finite element scheme for the Cahn-Hilliard equation ut + ∆(ε∆u − εf(u)) = 0 arising from phase transition in materials science, where ε is a small parameter known as an “interaction length”. The primary goal of this paper is to establish some useful a priori error estimates for the proposed numerical methods, in particular, by focusing on the dependence of the error bounds on ε. Quasi-optimal order error bounds are shown for the semi-discrete and fully discrete schemes under different constraints on the mesh size h and the local time step size km of the stretched time grid, and minimum regularity assumptions on the initial function u0 and domain Ω. In particular, all our error bounds depend on 1 ε only in some lower polynomial order for small ε. The cruxes of the analysis are to establish stability estimates for the discrete solutions, to use a spectrum estimate result of Alikakos and Fusco [3] and Chen [15], and to establish a discrete counterpart of it for a linearized Cahn-Hilliard operator to handle the nonlinear term on the stretched time grid. It is this polynomial dependency of the error bounds that paves the way for us to establish convergence of the numerical solution to the solution of the Hele-Shaw (Mullins-Sekerka) problem (as ε ց 0) in Part II [26] of the series.